Crack roughness and avalanche precursors in the random fuse model 



O 
O 



(N 



o 



CO 

o3 



i 

O 

o 



> 

00 

in 
i> 
o 

o 



-a 
c 

o 
o 



X 
S3 



Stefano Zapperi 
INFM UdR Roma 1 and SMC, Dipartimento di Fisica, 
Universita "La Sapienza" , P.le A. Moro 2, 00185 Roma, Italy 

Phani Kumar V.V. Nukala and Srdan Simunovic 
Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Rio 



je, TN 37831-6359, USA 



We analyze the scaling of the crack roughness and of avalanche precursors in the two dimensional 
random fuse model by numerical simulations, employing large system sizes and extensive sample 
averaging. We find that the crack roughness exhibits anomalous scaling, as recently observed in 
experiments. The roughness exponents (£, C,ioc) and the global width distributions are found to be 
universal with respect to the lattice geometry. Failure is preceded by avalanche precursors whose 
distribution follows a power law up to a cutoff size. While the characteristic avalanche size scales 
as so ~ L D , with a universal fractal dimension D , the distribution exponent r differs slightly for 
triangular and diamond lattices and, in both cases, it is larger than the mean-field (fiber bundle) 
value t — 5/2. 



I. INTRODUCTION 

Understanding the scaling properties of fracture in dis- 
ordered media represents an intriguing theoretical prob- 
lem with some technological implications ,1]. Experi- 
ments have shown that in several materials under differ- 
ent loading conditions the fracture surface is self-affmc 
and the out of plane roughness exponent displays a 
universal value irrespective of the material studied [3- 
In particular, experiments have been done in metals [4|, 
glass [fj| , rocks and ceramics Q , covering both ductile 
and brittle materials. 

It was later shown that the roughness exponent con- 
ventionally measured describes only the local properties, 
while the fracture surface instead exhibits anomalous 
scaling |10| : the global exponent describing the scaling 
of the crack width with the sample size is larger than 
the local exponent measured on a single sample It 
is thus necessary to define two roughness exponents a 
global (() and a local Q oc . Only the latter appears to be 
universal with a value Q oc ~ 0.8 Q. For the purpose of 
this paper, it is important to mention that experiments 
performed in quasi two-dimensional geometries, in wood 
[8j or paper yield a self-affine exponent close to the 
minimum energy surface result Q oc — 2/3. 

Scaling is also observed in acoustic emission experi- 
ments, where the distribution of pulses decays as a power 
law over several decades. Experimental observations have 
been reported for several materials such as wood [l2| . 
cellular glass concrete and paper ^f|, but uni- 
versality in the scaling exponents does not appear to be 
present. 

The experimental observation of scaling behavior sug- 
gests an interpretation in terms of critical phenomena, 
but a complete theoretical explanation has not been 
found. The motion of a crack front has been modeled as 
a deformablc line pushed by the external stress through a 
random toughness landscape. Deformation of the crack 
surface is caused by disorder and opposed by the elas- 



tic stresses. In certain conditions, the problem can be 
directly related to models and theories of interface de- 
pinning in random media and the roughness exponent 
computed by numerical simulations and renormalization 
group calculations [16L Il7| . Unfortunately, the numerical 
agreement between this theoretical approach and exper- 
iments is quite poor 0, 0, E(j ■ 

One aspect missing from the crack line model is the 
nucleation of voids in front of the main crack, an effect 
that has been shown to occur experimentally [2]]]. In 
this perspective, disordered lattice models appear to be 
more appropriate to describe the phenomenon. In these 
models the elastic medium is described by a network of 
springs with random failure thresholds. In the simplest 
approximation of a scalar displacement, one recovers the 
random fuse model (RFM) where a lattice of fuses with 
random threshold are subject to an increasing external 
voltage HI Him El El- The model has been numer- 
ically simulated to obtain the roughness of the fracture 
surface in two [53, El] and three dimensions E^E^I- The 
measured roughness exponents turn out to be similar to 
the one describing a minimum energy surface (or a di- 
rected polymer in d = 2) suggesting that crack formation 
occurs by an optimization process. 

In addition, the fracture of the RFM is preceded by 
avalanches of failure events [3l|, l.32l |33| . These are remi- 
niscent of the acoustic emission activity observed in ex- 
periments. The distribution of avalanche sizes (i.e. the 
number of bonds participating in an avalanche) follows a 
power law. In previous simulations the exponent resulted 
to be close to r = 5/2, the value expected in the fiber 
bundle model (FBM) [3ll35|. In the FBM, load is redis- 
tributed equally in all the fibers, representing thus a sort 
of mean-field limit of the RFM |33| ■ The load transfer in 
the RFM is long-ranged and is thus possible that RFM 
and FBM display universal behavior |36j| . An interme- 
diate case is provided by FBM with long-range (power 
law) load transfer : the difference with the RFM lies 
in the anisotropic current transfer function [36| . 
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Numerical simulation of fracture in the RFM is often 
hampered by the high computational cost associated with 
solving a new large set of linear equations every time a 
new lattice bond is broken. Previously, this fact has re- 
stricted the simulations to smaller lattice sizes and fewer 
statistical sampling of data, thereby affecting the quality 
of the results. Here, thanks to the new algorithm dis- 
cussed in Ref. 26], we report results of numerical simu- 
lations for large two-dimensional lattices (triangular and 
diamond) with extended statistics. We concentrate on 
the roughness of the final crack and the avalanche statis- 
tics preceding failure. 

Using local and global measurements for the roughness 
we find that cracks in the RFM follow anomalous scaling 
flfT ]. The local roughness exponent is found to be in the 
range C,i oc = 0.70 — 0.75 , while the global exponent falls 
in the range ( = 0.80 — 0.85. Although the difference 
between £ and Q oc is small it appears to be systematic. 
The results are obtained using the local width and the 
power spectrum methods and appear to be universal with 
respect to the lattice type. As a further test, we compute 
the width distribution that can be collapsed into a unique 
curve for different lattice sizes and types |38|. 

Next, we consider the distribution of avalanche sizes. 
The avalanche signal is not stationary and as the cur- 
rent is raised avalanches becomes larger and larger. The 
last avalanches, producing the failure of the sample, is 
typically much larger than the previous one and it fol- 
lows a normal distribution with a typical value scaling 
as s m ~ L 1A . Preceding avalanches are distributed as a 
power law with a cutoff increasing with the current. Inte- 
grating the distribution over all the values of the current, 
we find a power law up to a cutoff, scaling with the lattice 
size as L D , where D ~ 1.18 does not depend on the lat- 
tice type and is thus universal. The exponent describing 
the decay of the distribution is found instead to differ for 
triangular and diamond (square lattice with 45 degrees 
inclined bonds to the bus bars) lattices with a value which 
is always larger than the FBM value r = 5/2. 

The paper is organized as follows: in section II we de- 
fine the model, in section III we report the results on 
the crack roughness and section IV is devoted to the 
avalanche statistics and in section V we conclude. 



II. THE RANDOM FUSE MODEL 

In the RFM , the lattice is initially fully intact with 
bonds having the same conductance and random break- 
ing thresholds t, uniformly distributed between and 1. 
The burning of a fuse occurs irreversibly, whenever the 
electrical current in the fuse exceeds breaking threshold 
t of the fuse. Periodic boundary conditions are imposed 
in the horizontal direction to simulate an infinite system 
and a constant voltage difference, V, is applied between 
the top and the bottom of lattice system bus bars. Nu- 
merically, a unit voltage difference, V = 1, is set between 
the bus bars and the Kirchhoff equations are solved to 



determine the current flowing in each of the fuses. Sub- 
sequently, for each fuse j, the ratio between the current 
ij and the breaking threshold tj is evaluated, and the 
bond j c having the largest value, max 3 -p is irreversibly 
removed (burnt). The current is redistributed instanta- 
neously after a fuse is burnt implying that the current 
relaxation in the lattice system is much faster than the 
breaking of a fuse. Each time a fuse is burnt, it is nec- 
essary to re-calculate the current redistribution in the 
lattice to determine the subsequent breaking of a bond. 
The process of breaking of a bond, one at a time, is re- 
peated until the lattice system fails completely. At this 
point we analyze the morphology of the spanning crack. 

The same breaking sequence is obtained by raising the 
voltage difference or the total current at an infinitesimal 
rate. Doing this one can identify an avalanche as the 
set of fuses breaking between two successive increases of 
the voltage (or the current). In this paper, we follow 
Ref. [32] ] , considering only current driven avalanches. 
The avalanche size is defined as the number of fuses in 
an avalanche. 

Simulations are performed on two dimensional triangu- 
lar and diamond lattices of linear sizes going from L = 16 
up to L = 1024 (for the triangular lattice) or up to 
L = 256 (for the diamond lattice). The total number 
of bonds in the lattice is given by N = (31/ + 1)(L + 1) 
for the triangular lattice and N = 2L(L + 1) for the di- 
amond lattice. Several results discussed in the following 
sections could only be obtained under an extensive statis- 
tical sampling. Due to numerical limitations this could 
not be achieved for the largest lattice sizes. Each nu- 
merical simulation was performed on a single processor 
of Eagle (184 nodes with four 375 MHz Power3-II pro- 
cessors) supercomputer at the Oak Ridge National Labo- 
ratory. The statistically independent N con fi g number of 
configurations were simulated simultaneously on number 
of processors available for computation (the actual values 
of N con fig is reported in Table 1 of Ref. HJJ). 

III. CRACK ROUGHNESS 

After the sample has failed we identify the final crack, 
an example of which is reported in Fig. ^ The cracks 
typically display some limited amount of dangling ends 
and overhangs. We remove them and obtain a single 
valued crack line y x , where the values of x G [0, L] 
depend on the underlying lattice topology. Several 
methods have been devised to characterize the rough- 
ness of an interface and their reliability has been tested 
against synthetic data |4l|. If the interface is self-affine 
all the methods should yield the same result in the 
limit of large samples. For instance the local width, 

w ( l ) = iY^xiyx ~ (V0Ex») 2 ) 1/2 > where the sums 
are restricted to regions of length I and the average is 
over different realizations, should scale as w(l) ~ 1$ for 
I <C L and should saturate to a value W = w(L) ~ 
corresponding to the global width. The power spec- 



3 



trum S(k) = (y k y-k), where y k = ^2 x y x expi(2'Kxk/L), 
should decay as S(k) ~ k~^ +1 \ 

While numerical estimates with the two methods above 
could yield different results, it is also possible that the 
scaling is anomalous [Tof. This has been observed not 
only in various growth models 01 but also in fracture 
surfaces in granite and wood samples|llj. Anomalous 
scaling implies that the exponent describing the system 
size dependence of the surface differs from the local ex- 
ponent measured for a fixed system size L. In partic- 
ular, the local width scales as w(l) ~ .Z>~£ loc , so 
that the global roughness W scales as with Q > (i oc . 
Consequently, the power spectrum scales as S(k) ~ 

fe -(2Cioc+l)£2(C-Cioe). 

Previous measurements of the crack roughness in the 
two-dimensional random fuse model have been obtained 
studying the global roughness and anomalous roughness 
could not be detected. Here, thanks to the improved 
statistics and system size range, we reveal clear indica- 
tion of anomalous scaling behavior. In Fig. [21 we report 
the local width for diamond and triangular lattices for 
different sizes L. The curves for different system sizes are 
not overlapping even for I <C L as expected for anoma- 
lous scaling. The global width scales with an exponent 
C = 0.80 ± 0.02 and C = 0.83 ± 0.02 for diamond and 
triangular lattices respectively. On the other hand the 
local width increases with a smaller exponent, that can 
be estimated for the larger system sizes as (i oc ~ 0.7 for 
both lattices. A more precise value of the exponents is 
obtained from the power spectrum, which is expected to 
yield more precise estimates. Fig. [^reports the data col- 
lapse of the power spectra for different system sizes. The 
data are collapsed using C~C/oc = 0.1 and C — Ooc = 0.13 
for diamond and triangular lattices, respectively. A fit 
of the power law decay of the spectrum yields instead 
Cloc = 0.7 and Q oc = 0.74 for the two lattices, implying 
£ = 0.8 and Q — 0.87. The results are close to the real 
space estimates and we can attribute the differences to 
the bias associated to the methods employed |4l| . 

Although the value of £ — 0oc is small, it is significantly 
larger than zero so that we would conclude that anoma- 
lous scaling is present. While the local exponent is close 
to the directed polymer value ( = 2/3, the global value is 
much higher. In addition, the presence of anomalous scal- 
ing would invalidate universality between directed poly- 
mers and fracture: directed polymers should not display 
anomalous scaling. As for the question of universality of 
the random fuse model crack roughness exponents, the 
values measured above are quite close to each other and 
the differences could be due to size effects. In order to 
have a further confirmation of this, we have analyzed the 
distribution P(W) of the crack global width. This distri- 
bution has been measured for various interfaces in models 
and experiments and typically rescales as |38| 

P(W) = P(W/(W))/(W), (1) 

where (W) ~ Z> is the average global width. The crack 
width distribution has been measured for the random 



fuse model with limited statistical sampling. We show in 
Fig. 0] that the distributions can be collapsed well using 
Eq. n for diamond and triangular lattices. The plot in 
Fig. [S] shows that the collapsed distribution for the two 
lattices superimpose, which we consider as a further indi- 
cation of universality. Finally, the width distributions are 
well fit by a log-normal distribution as shown in Fig. [SJ 

IV. AVALANCHES 

The qualitative behavior of the avalanche statistics is 
well understood in FBM, which can be solved exactly rep- 
resenting a mean- field version of the RFM 132J . The FBM 
can be formulated as a parallel set of fuses, with ran- 
dom breaking threshold, under a constant applied cur- 
rent I. Thus each fuse carry the same current = I/n, 
where n is the number of intact fuses. The FBM has been 
solved exactly ant it is known that there is a critical value 
/ = I c at which the bundle fails through a macroscopic 
avalanche. For I < I c fuses burn in smaller avalanches, 
whose sizes are distributed as 

p(s,n = s-ih(-s/s*), (2) 

with 7 = 3/2, and h(x) is a cutoff function. The cutoff 
size s* increases with the current and close to I c diverges 
as s* ~ (I c —I) 1 / a with a — 1. One can then integrate the 
distribution over all the values of the current, obtaining 
an P(s) ~ s~ T with r = 7 + a = 5/2. 

Here we study the statistical properties of the 
avalanches in the RFM. We can use the scaling laws es- 
tablished for the FBM as a reference, with additional 
complications due to finite size effects. In Fig. [S]we re- 
port the integrated avalanche distribution obtained for 
different lattice sizes. We observe a power law decay cul- 
minating with a peak at large avalanche sizes. As in 
the FBM, the peak is due to the last catastrophic event 
which can thus be considered as an outlier and analyzed 
separately. When the last avalanche is removed from the 
distribution the peak disappears (see Fig. |BJ). 

The avalanche size distribution, once the last event is 
excluded, is a power law followed by an exponential cutoff 
at large avalanche sizes. The cutoff size so is increasing 
with the lattice size, so that we can describe the distri- 
bution by a scaling form 

P(s,L)=s-Tg(s/L D ), (3) 

where D represents the fractal dimension of the 
avalanches. To take into account the different lattice ge- 
ometries, it is convenient to express scaling plots in terms 
of N rather than L 

P{s 1 N) = s- T g{s/N D/2 ). (4) 

A powerful method to test these scaling laws, extract- 
ing r and D, is provided by the moment analysis |39j |. We 
compute the qth moment of the distribution M q = (s q ) 
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and plot it as a function of N. This defines an exponent 
dq as M q ~ N a ". If the data follow Eq.Hthen a q = for 
q < t — 1 and <7 9 = D(q + 1 — r)/2 for q > r — 1. In or- 
der to measure cr g , we consider lattice sizes from L = 16 
to L = 128 since the statistical sampling for larger sizes 
is not adequate to estimate correctly the cutoff so- The 
data displayed in Fig. [7] show that indeed u q is linear 
in q at large q and vanishes for small q. The curves for 
triangular and diamond lattice do not coincide: the two 
lines are parallel, indicating that D is the same, but the 
intersection with the x axis differs. By a linear fit we 
obtain t = 2.75 and D/2 = 0.59 for diamond lattices 
and t = 3.05 and D/2 = 0.585 for triangular lattices. 
To confirm these results we perform a data collapse us- 
ing the estimated values of the exponents and result is 
reported in Fig. [S] While the data collapse for diamond 
lattice is nearly perfect, some deviations are noticeable 
for the triangular lattice. 

From the analysis discussed above, we would conclude 
that the avalanche fractal D dimension is universal, but a 
significant difference is present for the exponent r. This 
difference could be due to lattice finite size effect as we 
will discuss later. In addition, the value of r appears to 
be larger than the mean- field result r = 5/2 obtained in 
the FBM. On the basis of less accurate results, it was 
conjectured in Ref. [SsJ that avalanches in the random 
fuse model are ruled by mean-field theory. The present 
results seem to rule out this possibility. 

So far we have considered avalanche statistics integrat- 
ing the distribution over all the values of the current. We 
have noticed, however, that the avalanche signal is not 
stationary: as the current increases so does the avalanche 
size. In particular, the last avalanche is much larger than 
the others. Its typical size grows as s m ~ N b , with 
b ~ 0.7 see Fig. 4 of Ref. |4jj ( s m is referred as rif — n p 
in that paper), while the distribution is approximately 
Gaussian as shown from the data collapse reported in 
Fig. H 

In Fig. EJ] we report the distribution of avalanche sizes 
sampled at different values of the current /. For each 
sample, we normalize the current by its peak value I c 
and divide the I* = I/I c axis into 20 bins. We then 
compute the avalanche size distribution p(s,I*) for each 
bin and average over different realizations of the disor- 
der. In Fig. 1101 we report this distribution for a diamond 
lattice of size L = 128. The distribution follows a law of 
the type 

p(s,n = ^ 7 exp(-s/s*), (5) 

with 7 ~ 1.9, while in the FBM 7 = 3/2. 

In order to extract the dependence of the cutoff s* 
on I*, we compute the second moment of the distri- 
bution (s 2 ). According to Eq. this should scale as 
(s 2 ) = (s*) 3 ~ 7 . Assuming that for large systems s* ~ 
(1 - I*)- 1 /" 7 (in the FBM this holds with a = 1), we 
expect that the singularity is rounded at small L as 

L D 

s * ~ (i-i*y/°L D + c' (6) 



where C is a constant. The second moment can be col- 
lapsed very well under this finite size scaling assumption 
with l/cr = 1.4 and D = 1.18 as shown in Fig. 1111 for 
the diamond lattice. The data collapse is consistent with 
the finite size scaling of the integrated distribution with 
a cutoff increasing as so ~ L D . In fact integrating Eq. 
we obtain 

P(s, L) ~ s-^+^ exp[sC/L D }. (7) 

which implies r = 7 + a. Using the estimated data we 
would obtain 7 + a ~ 2.6 in reasonable agreement with 
the integrated distribution result r = 2.75. 

We have performed the same analysis for the triangu- 
lar lattice, where we find similar scaling laws with 7 ~ 2 
and a = 1.3. This would give t = 2.7 that is quite off 
from the integrated distribution result r = 3.05. These 
variations could indicate some systematic error present 
in the triangular lattice results. We notice that while in 
the diamond lattice, at the beginning of the simulation, 
all fuses carry the same current, in the triangular lattice 
only two thirds of the fuses carry a current. As fuses 
break the current is redistributed becoming inhomoge- 
neous so that at breakdown this lattice effect should not 
be visible. In fact scaling exponents computed at failure, 
like the roughness exponent or the avalanche cutoff, do 
not depend on the lattice type. On the other hand, the 
integrated avalanche distribution is affected by the entire 
rupture process and the estimated exponent could thus 
be biased. 



V. CONCLUSIONS 

In this paper we have revised some statistical prop- 
erties of fracture in the random fuse model using an 
improved statistical sampling and larger lattices than 
what previously done in the past. We have analyzed 
the roughness of the final crack for diamond and trian- 
gular lattices. The local roughness exponent is found 
to be Q oc — 0.72 ± 0.03 and appears to be different 
from the global roughness exponent which turns out to 
be C = 0.83 ± 0.04. These results have been obtained 
from the local width and the power spectrum methods 
and the error bars above merely represent the spread 
of the estimated exponents using various methods and 
lattice types. The data suggest that anomalous scal- 
ing is present, as already found in fracture experiments 
. The numerical value of the local exponent is in rea- 
sonable agreement with the experiments on quasi two- 
dimensional materials 8, 9]. As a further test for uni- 
versality, we have also evaluated the width distribution 
|S8j that can be collapsed into a single curve for differ- 
ent lattice sizes and types. From the theoretical point of 
view, our results seem to exclude the minimum energy 
surface exponent of ( = 2/3. While the local exponent 
is close to that value, the global exponent is definitely 
higher. In addition anomalous scaling is not expected 
for that model. Thus the origin of measured roughness 
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exponents and its theoretical explanation remains still 
open. 

We have also analyzed the scaling of failure precur- 
sors, computing the distribution of avalanche sizes. The 
extensive statistical sampling employed allowed us to ob- 
serve a power law decay up to a cutoff, which was not 
visible in previous simulations 

HEIl^l The cutoff size is 
found to increase with the lattice size as so ~ L D , where 
the exponent D ~ 1.18 depends very little on the lattice 
size. It is interesting to notice that for self-affine lines of 
roughness, £, one expects a fractal dimension D = 2 — £ 
|42|. If we plug into this expression the global rough- 
ness results obtained above for the final crack, we obtain 
D ~ 1.13 — 1.20. This could imply that the geometri- 
cal properties of the precursors are the same as that of 
the final crack. On the other hand, the exponent of the 
avalanche size distribution displays significant variations 
with the lattice type (i.e. r = 2.75 and r = 3.05 for 
diamond and triangular lattices respectively) and is sig- 
nificantly different from the mean- field result r = 5/2 
that was conjectured to be valid in j^] . 

The integrated avalanche distribution is due to the con- 
volution of the avalanche distribution measured at differ- 
ent values of the current. We have shown that the non- 
integrated distribution is given by a power law with an 



exponential cutoff that increases with the current. The 
combined analysis of the distribution with respect to cur- 
rent and lattice size can be performed using finite size 
scaling. The behavior of the model is similar to the FBM, 
as noticed in Ref. |32|, but the numerical values of the 
exponents change. For the diamond lattice we estimate 
7 = 1.9 and a = 1.4, while the FBM yields 7 = 3/2 
and (7 = 1. Similar results hold for the triangular lattice 
although the scaling there appears to be less clear. 

It would be interesting to understand these results the- 
oretically by the renormalization group, using the mean- 
field theory as a reference. Steps in this direction have 
been made in Ref. [3(| but the complicated (dipolar) 
structure of the current redistribution function makes 
the problem very hard to deal with. Long-range inter- 
actions appear to be crucial in the appearence of scaling 
behavior, since local fracture models yield abrupt failure 
without large precursors 43]. A similar scenario is char- 
acteristic of first-order phase transitions occuring close to 
a spinodal. In that case spinodal scaling is only seen in 
mean- field or with long range interaction |44j |. The fact 
that the exponents deviate from mean-field ones, how- 
ever, calls for a more detailed understanding of the origin 
of scaling in the random fuse model. 
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FIG. 1: The final crack in a triangular lattice of size L = 1024 (a detail is shown in the inset). The crack displays some dangling 
ends and overhangs that are removed before performing the analysis. 
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FIG. 2: The local width w(l) of the crack for different lattice sizes in log-log scale. A line with the local exponent Q oc = 0.7 is 
plotted for reference. The global width displays an exponent £ > Q oc . Data are shown for diamond (left) and triangular (right) 
lattices. 




FIG. 3: The power spectrum of the crack S(k, L) for different lattice sizes in log-log scale. The slope defines the local exponent 
as — (2£; oc + 1). The spectra for all of the different lattice sizes can be collapsed indicating anomalous scaling. Data are shown 
for diamond (left) and triangular (right) lattices. 
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FIG. 4: The distribution of crack width for different lattice sizes can be collapsed using their average value. Data are shown 
for diamond (left) and triangular (right) lattices. 
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FIG. 5: The distribution of crack width is universal for diamond and triangular lattices since all the curves can be collapsed 
together. A fit with a lognormal distribution is shown by a dashed line. 
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FIG. 6: The distribution of avalanche sizes for triangular lattices of different sizes. The peak at large size is due to the last 
avalanche, corresponding to catastrophic failure (right). On the left figure we show the same distribution without the last event 
and with logarithmic bins. 
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FIG. 7: The exponent a q ruling the scaling of the gth moment for triangular and diamond lattice. The shift in the lines 
indicates a difference in the value of r. 
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FIG. 9: Data collapse of the distribution of the last avalanche for diamond (left) and triangular (right) lattice. 
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FIG. 10: The avalanche size distributions sampled over a small bin of the reduced current I* for a diamond lattice of size 
L = 128. The dashed line represents a fit according to Eq.^with 7 = 1.9. 
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FIG. 11: The second moment of the avalanche size distribution as a function of the reduced current 1 — I* for diamond lattices 
of different sizes (inset). The curves can be collapsed using the finite size scaling assumption reported in Eq. [SJwith 7 = 1.9, 
D = 1.18 and 1/a = 1.4. 



